Postavljanje jednacina
syms ug uC uL L C DiL DuC m u iC iL Um t
jednacine = [ug == uC + uL,...
iC==iL]
jednacine = 
Eliminacija i redukcija sistema
sistem = eliminate(jednacine, [uL,iC,u])
sistem = 
Resavanje po izvodima stanja
resenje = solve(sistem, [DiL,DuC])
resenje =
DiL: -(uC - ug)/L
DuC: iL/C
Resavanje sistema diferencijalnih jednacina
jednacineStanjaFun = subs([diff(uC)==resenje.DuC; diff(iL)==resenje.DiL], {iL,uC}, {iL(t),uC(t)})
jednacineStanjaFun(t) =

zamene = ug==Um*sin(t/sqrt(C*L))*heaviside(t)
zamene =

jednacineStanjaFunR=subs(jednacineStanjaFun,lhs(zamene),rhs(zamene))
jednacineStanjaFunR(t) =

resenjeDiff=dsolve(jednacineStanjaFunR,[uC(0)==0,iL(0)==0])
simplify(resenjeDiff.iL)
ans =

numzamene=[L==1e-4, C==1e-6, Um==1, m==1]
numzamene =

fplot(t, subs(resenjeDiff.iL, lhs(numzamene), rhs(numzamene)),[0 0.002])